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We examine the time-dependent behavior of a nonlinear system driven by a two-frequency forcing. 

By using a non-perturbative approach, we are able to derive an asymptotic expression, valid in the 
long-time limit, for the time average of the output variable which describes the response of the 
system. We identify several universal features of the asymptotic response of the system, which are 
independent of the details of the model. In particular, we determine an asymptotic expression for 
the width of the resonance observed by keeping one frequency fixed, and varying the other one. We 
show that this width is smaller than the usually assumed Fourier width by a factor determined by 
the two driving frequencies, and independent of the model system parameters. Additional general 
features can also be identified depending on the specific symmetry properties of the system. Our 
results find direct application in the study of sub-Fourier signal processing with nonlinear systems. 
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I. INTRODUCTION 

Nonlinear systems driven by time-dependent forcings 
give rise to a wealth of intriguing phenomena, from chaos 
to the stochastic resonance phenomenon and the ratchet 
effect [T]. Given the complexity of the systems, a com¬ 
plete analysis of these phenomena is usually only possible 
via a numerical treatment, which is inherently model de¬ 
pendent. A more general analysis, less model dependent, 
is usually possible only for the stationary state of the 
system, in the presence of system symmetries. These ap¬ 
proaches thus do not easily allow the determination of 
the asymptotic scaling properties of the systems, whose 
identification is essential for the understanding of unique 
features of non-linear systems, as the recently highlighted 
sub-Fourier dynamics which may open the possibility of 
sub-Fourier information processing Hi- 

In this work, we study theoretically the time- 
dependent behavior of a generic nonlinear system driven 
by a two-frequency forcing. To specify the notation, we 
refer to a single-particle model, and assume that the av¬ 
erage particle velocity is the quantity of interest for the 
study of the response of the system. This is for exam¬ 
ple the case relevant to the study of directed motion in 
driven systems. By using a nonperturbative approach, 
we are able to derive an asymptotic expression, valid in 
the long-time limit, for the average particle velocity. This 
expression allows us to identify several universal features 
of the asymptotic response of the system, which are inde¬ 
pendent of the details of the model. Additional general 
features can also be identified depending on the specific 
symmetry properties of the system. 


This paper is organized as follows. In Sec. [TTj we in¬ 
troduce the system of interest and state the problem un¬ 
der consideration. In Sec. m we establish some general 
asymptotic properties of the nonlinear system, valid in 
the long-time limit. In Sec. IV we illustrate our theo¬ 
retical results with the help of a one-dimensional model 
consisting of a Brownian particle, moving in a periodic 
potential, under the influence of a biharmonic force. We 
also investigate how the specific symmetries of the evo¬ 
lution equation affect our asymptotic results. Finally, in 
Sec. |Vj we present conclusions for the main findings of 
our work. 


II. STATEMENT OF THE PROBLEM 

Let us consider a generic single-particle model in which 
the time evolution of the particle’s position vector r{t) 
is described by a differential equation, which may be 
stochastic or deterministic depending on the presence or 
absence of noise sources. From the solution of this dif¬ 
ferential equation with an appropriate set of initial con¬ 
ditions at an arbitrary time to, one can determine the 
particle’s position vector at any later time tf > to- The 
average velocity between these two times is then defined 
as 




t{ 

2 y dt (r(t)) 
*0 


(r(tf) - r(to)) 
Te ^ 


( 1 ) 


where Tg = t{ — to is the elapsed time, the overdot rep¬ 
resents time derivative, and (...) the ensemble average 
with respect to the noise, if present. The infinite-time 
average velocity is obtained by taking the limit 


V = lim v*o,to+T„ ^ 
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In most dissipative systems this limit exists and is inde¬ 
pendent of the initial time to and of the initial conditions. 

In simulations, or real experiments, the variable r(t) is 
known only in a finite time-interval. Thus, the limit in 
Eq. § can be calculated only approximately by choos¬ 
ing a sufficiently large value for the elapsed time Tg. To 
be more precise, Tg must be chosen to be much larger 
than all other relevant time-scales in the problem. In 
particular, if the system is subjected to a periodic time- 
dependent forcing function with period T, then a neces¬ 
sary condition is that Tg ^ T. 

The main focus of the present work is to determine 
scaling properties with the finite-time Tg. Our results will 
be applicable to any system, with a well-defined infinite¬ 
time average § , containing time-dependent forcing func¬ 
tions of the form 

F{t) = ^ . ( 3 ) 

In the above expression, $ {6i, O2) is a vector function [ 3 ] 
periodic with respect to both arguments with period 
27 r [i.e., $( 6 >i-k 27 r, 6 » 2 ) = $ (0i, 6»2-f 27 r) = $(6li,6>2) 
^(^11^2) G and uJi and UJ2 are two positive parame¬ 
ters with the dimensions of frequency. Bifrequency forc¬ 
ing functions of this type have been extensively used in 
the literature (e.g., [ 310 [ 6 ]). If wi and 102 are commen¬ 
surate, i.e., LO2 = poJi/q with p and q coprimes, F(t) is 
periodic of period T = 2 T:q/uJi. 

To be more specific, we consider the following sit¬ 
uation: one of the frequencies, say wi, is kept fixed, 
while the other frequency, W2, is scanned around a fixed 
value, a;2,0) which is commensurate with wi. This corre¬ 
sponds to quantum-optical setups of recent experiments 
[ 3 ]. We will indicate the ratio between 102,0 and loi with 
^2,0/^1 = Po/qoj with po and qo coprimes. We aim 
to derive the general asymptotic properties of the func¬ 
tion v‘“’‘‘’+^=(a;2) in an interval around uj2,o- As a first 
step, we analyze to which extent the infinite-time average 
V(a;2) can be approximated by the finite-time average 
v*°’‘°+^'’(w2). This analysis will reveal a general scaling 
property which will allow us to derive the asymptotic 
limit. 

As discussed above, a necessary condition to approx¬ 
imate the infinite-time limit V(w2.o) by the finite-time 
average v*“>‘“+^®(a;2,o) is that Tg ^ Tq = 2 -KqoluJi] i.e., 
the elapsed time should be significantly larger than the 
period of the driving. However, as we now show, no 
matter how large Tg is chosen, there always exists a fre¬ 
quency interval centered at a;2,o within which it is not ap¬ 
propriate to approximate the average velocity V(a;2) by 
the finite-time average except for the value 

UJ2 = 102,0 itself. Indeed, let us consider the frequency in¬ 
terval defined by the inequality \i02 — W2,o| < 27 r/(go 7 e), 
and choose any value a;2,i ^ 102,0 in this interval of 
the form a;2,i = piioi/qi, with pi and qi being co¬ 
prime integers. Taking into account that there are in¬ 
finitely many rational numbers within any real inter¬ 
val of nonzero length, the number of possible choices 
for 102,1 is infinity. Since \piqo — Poqi\ > Ij it fol¬ 


lows that ioi/{qoqi) < \oo2,i — a;2,o| < 27 r/(goFe). Thus, 
Ti = 2t: qijioi > Tg and, consequently, for the chosen 
value of Tg, it is not appropriate to approximate V(a;2p) 
by v*°’*°+^'’(a;2,i). The same conclusion can be drawn 
if i02,\loo\ is irrational, since any irrational number can 
be approximated to any degree of accuracy by rational 
numbers. 

The above result shows that, given an averaging time 
Tg, there is a frequency interval of length 47r/((7o?e) cen¬ 
tered at i02.o in which the behavior of the infinite-time 
limit V(a;2) cannot be inferred from that of the finite¬ 
time average v*°’‘°+^'(a;2). It is precisely the scaling 
with Tg of the length of such an interval which allows 
us to derive the asymptotic behavior of the function 
(a;2), the quantity that is directly measurable in 
a simulation or experiment, within the above-mentioned 
range of 102 values. 


III. DERIVATION OF THE ASYMPTOTIC 
PROPERTIES 

In this section, we analyze the asymptotic behav¬ 
ior of the average velocity for 102 values 

such that \uj2 — (02,o\ is of the same order of magnitude 
as 2 'K/{qoTg), and sufficiently large values of Tg. To 
this purpose, we make a change of variables from the 
variable 8102 = i02 — (02,0 to a dimensionless variable 
A(02, which is expressed in terms of 6102 according to 
A(02 = qoS(02Tg/( 2 tt). The average velocity as a function 
of this new variable will be given by (AW2) = 

Y*o,to+T ^022 0 + 2'KACo2/{qoTe)], where the subscript a;2,o 
indicates the midpoint of the region under study. In order 
to enforce the correct scaling derived in the previous sec¬ 
tion, we will consider the limit as Te —>■ 00 and 8102 —>■ 0, 
so that Aw 2 is held constant. Then, the leading-order of 
the aforesaid asymptotic behavior can be obtained from 
the limit 

(AA2) = ^ liin (AA2), ( 4 ) 

ie —>-+00 ’ 

while keeping Aw2 fixed. As a fact useful in the following, 
we notice that the time-dependent forcing function in 
Eq. ([^ can be rewritten in the form 

F{t) = $ [lOit, I02,ot + (p{t)] , ( 5 ) 

where ip{t) = 8i02t = 2t:AC b2t/{qoTg). 

In order to find a more explicit expression for the 
asymptotic result Vui2,oi^'^2), it is convenient to intro¬ 
duce an auxiliary function of Tg, N(Tg), with the follow¬ 
ing three properties [ 7 ]: (i) N{Tg) is dimensionless, inte¬ 
ger, and positive for all Tg > 0, (ii) limT,->.+oo N{Tg) = 
-boo, and (iii) liuiT^^+oo N{Tg)/Tg = 0 . By using this 
function, one can divide the interval into N{Tg) 

subintervals of equal length Ar(Te) = Tg/N{Tg). Then, 
taking into account the definition in Eq. Q, it is possible 
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to rewrite Eq. Q as 


_/V(Te) 

V.,.„(Aw2)=^lim ^ (6) 

ie —^+00 1 q ' ^ ’ 

i=i 

where r^- = Tj{Te) = to + jAT{Te). 

We now use two properties that follow from the defini¬ 
tion of the function N{Te). First, from (hi) it follows that 
lim7’^_,._i_oo AT(Te) = -1-00, i.e., the length of the above- 
mentioned subintervals tends to infinity as Tg —>■ -|-oo. 
Second, from (ii) it follows that lim7’^_,.+oo Aip{Te) = 0 , 
where AipiTg) = 2 ttAC j2AT{Te)/{qoT^) is the increment 
in Lp{t) corresponding to an increment Ar(Te) in t. In 
the calculation of this last limit, the variable Aui2 is held 
constant. Thus, within each of these subintervals the 
function ip{t) in Eq. ([^ can be approximated by a con¬ 
stant (for instance, by its value at the left endpoint of 
the subinterval). This last approximation becomes exact 
in the limit Te —>■ -l-oo. From these two properties and 
Eq. 1 ^, it follows that 

Af(Te) 

^ V(a;2.o,<p,-i)A^, 

i=i 

( 7 ) 

where Aip = Aip{Te), ipj = <fj{Te) = (p[Tj{Te)] = 
27 rAw 2 <o/('?o 7 "e) + jAip(Te), and 'V{ui2,o,(pj-i) is the 
infinite-time average velocity corresponding to a forcing 
function of the form ([^ but with a time-independent 
phase shift equal to Finally, taking into account 

that the sum in Eq. Q is a Riemann sum, as well as, 
that limT,-).-i-oo <^0(71!) = 0 and limr^_>+oo (/?Ar(Te) ('7'e) = 
2TrAu)2/<lo, one obtains 

2'kAu}2IQ0 

V^^ „{ACj2) = 2 JAQ 2 I • ( 8 ) 

0 

This expression is one of the central results of the present 
work. It shows that the asymptotic behavior V^2,o(^^2) 
is completely determined by 'V{uj2,o, ‘p), that is, by the 
infinite-time average velocity corresponding to a forcing 
function Fo(t, (/?) = ^ (wit, a;2,ot + (/s), considered as a 
function of the time-independent phase shift p. 

Equation ^ allows us to draw conclusions about the 
behavior of the asymptotic result Vuj2,oA^‘i) from the 
specific features of the function V(a;2,0j p)- For instance, 
it is easy to show that, under quite general assumptions, 
V(a;2,o, p) is a periodic function in (p with period 27 r/go; 
i.e.. 


V(a;2,o,V^) = V(a;2.o,<P + 27r/qo) (9) 

for all if. Indeed, since po and qo are coprime, from the 
Bezout identity [8] it follows that there exist integers I 
and n such that Ipo+nqo = 1 . From this equality, one can 
readily check that Fo{t + 2 nl/oJi,(p) = Fo{t,(p + 2 Tr/qo). 
Therefore, by assuming that all the time-dependent 


noise sources (if any) can be described by stationary 
stochastic processes, one has that v‘o’*o+'^=(a;2,o, ^z?) = 
v*°’‘‘’+^'(a;2,0)'^ _l_ 27 r/qo), with = to + 27 rl/a;i. Equa¬ 
tion ([^ then follows from the independence of the limit 
in Eq. ([^ on the initial conditions. Notice that, accord¬ 
ing to Eq. the average of 'V{u}2,o,p) over a period 
27 r/( 7 o is sirn^y given by Vu2oi^)- Taking into account 
Eqs. (|^ and it is easy to verify that 

Vc...o(fc) = V....o(l) (10) 

for any nonzero integer k. Similarly, one can check that 


{-k/2) - V.,,„(l) = V^,,„(l) - {k/2) (11) 


for any odd integer k. Equations ( 10 ) and 0 re¬ 
veal two remarkable properties of the asymptotic result 
Vcj^ p(Aw 2) which follow directly from the periodicity of 
the infinite-time average velocity expressed in Eq. (§. 

In order to gain some insight into the qualitative fea¬ 
tures of the function p(Aw2), let us consider first its 
behavior in the vicinity of the point Aw2 = 0 . From 
Eqs. (j^ and (^, it is clear that Vtj2,o(9) = V(w2,o,0). 
Furthermore, by appljdng L’Hopital’s rule, one can eas¬ 
ily show from Eq. ® that limAia-s-o V(.,j2_o (^‘7’2) = 
V(a;2,Oj0). Thus, Vi^2g(AL02) is continuous at Aw2 = 0 
and approaches the value V(w2,0) 0) as Acj2 tends to zero. 
To understand what happens when |A(I)2| increases, it 
is convenient to introduce the function 'D^^^^{Au}2) = 
27 rA£i 2 [V<^2,o(^‘ 7’2) - V^2.o(l)] / 9 o, so that 


V^2,o(^‘ 7>2) — VaJ 2 ,o(l) + 


go 1^012,oA'^2) 

2ttAu}2 


( 12 ) 


From the definition of T>^^,^{ACj2) and Eq. (|^, 
it is straightforward to see that T>,^^g{Auj2) = 

j2nAuj2jqo j-y ^ ( 1 )]. Therefore, given 

that V(a;2,0)'z2) is a periodic function in tp with pe¬ 
riod 2'K/qo and average value Vw2,o(l)’ il follows that 
'D^^g{Auj2) is a periodic function of Ad;2 with period 


1 . Consequently, according to Eq. ( 12 ), the asymptotic 
result V^2 o(^^ 2) approaches the value Vtj2o(l) ™ 
oscillatory manner as |Aa)2| increases, i.e., as UJ2 moves 
away from the value 012,0- The result given in Eq. ( 10 ) is 
a further indication of this oscillatory behavior. 


IV. MODEL-SPECIFIC EXAMPLES 

To illustrate our results, we will use a one-dimensional 
model consisting of a Brownian particle, with mass m 
and position a;, moving in a periodic potential U(x), with 
period A, under the influence of a biharmonic force of 
the form F(t) = A[cos(a;it) -I-cos(a;2t)]. The dynamics 
of this system is described by the stochastic differential 
equation 

mx(t) = — ax(t) — U' [a:(t)] -I- F(t) + ^(t ), ( 13 ) 










4 


where U'{x) is the derivative with respect to x oiU(x), a 
is the friction coefhcient, ^(t) is a Gaussian white noise of 
zero mean and autocorrelation (^(t)^(s)} = S{t — s), and 
r is the noise strength. The periodic potential U{x) can 
always be expressed as U{x) = Uo'^{2ttx/X)/2, where 
Uq is a constant with the dimensions of energy and 4 ^( 0 ) 
is a dimensionless periodic function with period 27 r. By 
expressing Eq. ( 131 in terms of the dimensionless quan¬ 


tities i = ujit and x = nx/X, the model is completely 
determined by five dimensionless parameters, namely, 
TO = mujifa, Uq = {X^aoJi), A = ttA!{X auji), 

r = TT^r /{X^a^uji), and Cj2 = 

In Fig. [2 we depict the dependence of the dimen¬ 
sionless average velocity i/°^^'(Aw2)/no on the variable 
Aa;2, where we have introduced the typical velocity 
Vg = Awi/tt. Since the model is one-dimensional, the 
boldface vector notation has been dropped. For the sake 
of simplicity, in this first figure we have considered the 
overdamped deterministic case (to = F = 0 ). Further¬ 
more, in order to avoid the presence of symmetries which 
are discussed later, we have chosen a highly asymmetric 
periodic potential (see figure caption for details). The 
results shown in Fig. clearly demonstrate that, for the 
elapsed time values considered (i.e., for Te = and 

Tg = the average velocity has es¬ 

sentially reached its asymptotic limit V^2 j,(Aa;2). Con¬ 
sequently, D°’J'^{Auj2) displays the typical oscillatory be¬ 


havior predicted by Eq. ( 12 ), and satisfies the condi¬ 
tions imposed by Eqs. (10) and ([Tl] ). It is to be noted 
that, in this figure, V,^ 2 ,o( 0 ) = ^(^o, 0 ) = A/Tg, where 
Tg = 27 rgg/a;i = Gtt/uji is the period of the biharmonic 
force (see Ref. [ 9 ] for an explanation of this fact). 

We return now to the general case to investigate 
how the specific symmetries of the evolution equation 
affect the shape of g(A£l 32 )- As a first example, 
let us assume that the differential equation describ¬ 
ing the evolution of r(t) in the presence of the forcing 
function Fg (t, ip) is invariant under the transformation 
Si : {r, r, r, t, ip} {ri - r, -r, -r, t + ti,ip- ir/qg}, 
where ri and ti are constant parameters. Then, if r(t) is 
a solution of this differential equation corresponding to 
a given value ip of the phase shift, R(t) = ri — r(t -|- ti) 
is also a solution corresponding to the displaced value 
ip + Tr/qg. Since (R.(t)) = —(r(t + G)) for all t, it follows 
from Eq. ([^ that V(w^,(p) = -V(a;2,c, -b Tr/gg) for 
all ip. Thus, from Eq. (j^, it is clear that V^., ( 1 ) = 0 . 
Consequently, in this case, Eqs. ( 10 ) and ( 11 ) can be 
rewritten, respectively, as Vi.c;2,o(^) ~ 0 fo'^ nonzero 
integer k, and V^ja.o (~fc/ 2 ) = ~ Vu;2,o (^/ 2 ) for any odd 
integer k. According to the first of these last two equa¬ 
tions it follows that, assuming that there exists a directed 
motion for a;2 = W2.0; this current can be suppressed 
by slightly shifting the frequency UI2 from the value a;2,o 
to the nearby values ClI2,o i 2 'K/{qgTg). This offers the 
possibility to discriminate two neighboring values of the 
frequency UJ2 differing in a quantity 2 Tr/{qgTg), which is 
qg times smaller than the usually assumed Fourier width 



FIG. 1. (Color online) Dependence of the dimensionless 
average velocity i'}}’^g{ALj 2 )/vo on the dimensionless vari¬ 
able Auj 2 for the model given by Eq. ( |13| ), with U{x) = 
Uo [cos(27ra;/A) -I- l/4cos(47ra:/A -I- tt/S)] /2. The value A 1 Z 12 = 
0 corresponds to 1 x 2 = 4)2 ,0 = Po/Qo, with po = 2 and 
qo = 3. The other parameters are ih = T = 0, Uo = A = 1, 
Te = 10'^ (circles) and Te = (red crosses). The 

solid line indicates the asymptotic result V,J 2 0 (A4’2)/'Co ob¬ 
tained by numerically evaluating the integral in Eq. ([^ . The 
dotted line shows the value Vcj 2 o(l)/'^o ~ —0.0633. The in¬ 
set depicts the dimensionless infinite-time average velocity 
V{uj2,o, ip) I Vo as a function of the phase shift ip. The infinite¬ 
time limit has been approximately calculated by evaluating 
the average velocity (uj 2 ,o, p) for Te = 


27T/Tg 0 . 

As a second example, let us assume now that the afore¬ 
said differential equation is invariant under the transfor¬ 
mation S2 : {r, r, r, t, ip} {r2 - r, f, -r, t2 - t, -ip}, 
with r2 and t2 being constants. Then, if r(t) is a so¬ 
lution corresponding to a given value ip of the phase 
shift, R(t) = r2 — r(t2 — t) is also a solution corre¬ 
sponding to the opposite value —ip. Thus, following 
a similar reasoning as above, it is easy to show that 
V(a;2m—1^) = V(a;2,c,(^) for all P- Consequently, from 
Eq. ^ it readily follows that V(.j2,o(^<1)2) is an even 
function of Aw2; i.e., (-AA2) = Vuj2,oi^^2) for all 

A£l 32 - In this case, Eq. ( 11 ) can be rewritten in the form 
VtcJ2,o (^/ 2 ) = Vcj2,o(l) for any odd integer k or, if the 
differential equation is also invariant under the transfor¬ 
mation Si, in the form Vui^ 0 ~ 

In Fig. we depict again the dependence of the di¬ 
mensionless average velocity i/°^^j(AA2)/ug on the vari¬ 
able Aui2 for the same model as in Fig. given by 
Eq. ( 13 ), but with a symmetric potential and a nonzero 
noise strength (see figure caption for details). We con¬ 
sider both the overdamped case (top panel) and the un¬ 
derdamped case (bottom panel). In both cases, it is easy 
to see that the corresponding evolution equations are in¬ 
variant under the one-dimensional version of the transfor¬ 
mation Si with, e.g., tl = —tt/wi and cci = 0 (TUI. Conse- 
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FIG. 2. (Color online) The same as in Fig. [^but for a sym¬ 
metric potential U{x) = Uocos{2'kx/\)/2, a dimensionless 
noise strength F = 0.05, and two values of the dimension¬ 
less mass, namely, rh = Q (top panel) and m = 1 (bot¬ 
tom panel). The rest of the parameters are the same as 
in Fig. The dashed lines indicate the results obtained 


from Eqs. ( |14[ ) (insets) and (151 (main figures). The model- 
dependent parameters are x ~ V{ijJ 2 ,o, x)ho = 0.073 

(top panel), and x ~ 0.963 and V{ijJ2fi,x)ho ~ 0.067 (bot¬ 
tom panel). These dashed lines are indistinguishable from the 
corresponding solid lines. The dotted lines show the value 
V<,;2,o(l)/«o = 0. 


quently, in both panels Vt.j2,o(l) = 0 (see the dotted lines). 
However, only in the case of the overdamped dynamics, 
there exist values for t2 and X2 such that the evolution 
equation is invariant under the one-dimensional version 
of the transformation 5 * 2 , e.g., t2 = 0 and X2 = A/ 2 . 
Thus, only the top panel shows the features associated 
with this transformation. 

At this point, it is important to remember that, in 
many one-dimensional systems driven by biharmonic 
forces, V{uj2^o,^p) can be expressed approximately as 

V (w2,o, <p) = V (w2.o, x) cos [go(<^ - x)] , ( 14 ) 


where y is the value of (p in the interval [ 0 , 27 r/qo) for 
which V(uj2,o,‘P) attains its maximum value V{uj2.o,x)- 
As shown in Ref. [IB, the functional form in Eq. ( | 14 [ ) 
is a consequence solely of the system symmetries, being 
independent of the interaction details. By introducing 
Eq. (jldl) in Eq. ( 1 ^, one obtains 


= V{uj2,0,x) 


sin ( 27 rAw 2 - Qox) + sin(gox) 


27rA£l’2 


(15) 


This ex pression is, in fact, a special case of 
Eq. ( 12 ) where V^2,o(l) = 0 and (,(A£l)2) = 
V (w2.o, X) [sin (27rAw2 - qox) + sin(gox)] /Qo- 

The significance of Eq. ( 15 ) lies in the fact that 
it is an explicit expression for the asymptotic result 
V^2 p(A£l’2) with only two model-dependent parameters, 
namely, y and V{uj2,o,x)- Particular cases of the sub- 
Eourier scaling with 27 r/(goTe) evidenced by Eq. ( 15 ) 
have been proposed in the literature based on heuristic 
arguments mm, as well as observed in cold-atom exper¬ 
iments [ 3 ]. The results obtained from Eqs. ( 14 ) and ( 15 ), 
with an appropriate choice of the model parameters, are 
indicated in Eig.j^with dashed lines. The agreement is so 
good that it is nearly impossible to distinguish between 
the solid and the dashed lines. By contrast, a glance at 
the inset of Fig. reveals that it is not reasonable to ap¬ 
proximate that curve by a function of the form given by 
Eq. (HtI). 


V. CONCLUSIONS 


In conclusion, we determined a universal asymptotic 
behavior of the average of a variable over a finite time in¬ 
terval. Our main results, expressed by Eqs. (|^-( 12 ), are 
applicable to any system driven by a bifrequency forc¬ 
ing function and with a well-defined infinite-time aver¬ 
age, independent of the initial conditions. With the help 
of Eq. (|^ and some symmetry considerations, we have 
obtained some qualitative and quantitative information 
about the aforesaid asymptotic behavior. In particular, 
we have derived an asymptotic expression for the width of 
the resonance observed by keeping one frequency fixed, 
and varying the other one. We have showed that this 
width is smaller than the one predicted by the Fourier 
inequality [ 2 ] by a factor determined by the two driving 
frequencies, and independent of the model system pa¬ 
rameters. The nonperturbative approach presented here 
to determine the asymptotic scaling finds direct applica¬ 
tion in the identification of non-linear systems displaying 
sub-Fourier resonances mia. 
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